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Luminescence in Ce doped materials corresponds to a transition from an excited state where the 
lowest Ce 5 c? level is filled (often called the (Ce"^^)* state) to the ground state where a single 4/ 
level is filled. We have performed theoretical calculations based on Density Functional Theory to 
calculate the ground state band structure of Ce-doped materials as well as the (Ce"^^)* excited state. 
The excited state calculations used a constrained occupancy approach by setting the occupation 
of the Ce 4/ states to zero and allowing the first excited state above them to be filled. These 
calculations were performed on a set of Ce doped materials that are known from experiment to be 
scintillators or non-scintillators to relate theoretically calculable parameters to measured scintillator 
performance. From these studies we developed a set of criteria based on calculated parameters that 
are necessary characteristics for bright Ce activated scintillators. Applying these criteria to about 
a hundred new materials we developed a list of candidate materials for new bright Ce activated 
scintillators. After synthesis in powder form one of these new materials (Ba2YCl7:Ce) was found to 
be a bright scintillator. This approach, involving first-principles calculations of modest computing 
requirements was designed as a systematic, high-throughput method to aid in the discovery of new 
bright scintillator materials by prioritization and down-selection on the large number of potential 
new materials. 

PACS numbers: 71.15.Qe, 71.20.Ps, 78.70.Ps 



I. INTRODUCTION 

Inorganic scintillators are extensively employed as ra- 
diation detector materials in many fields of applied and 
fundamental research such as medical imaging, high en- 
ergy physics, oil exploration, astrophysics and nuclear 
materials detection for homeland security and other 
applications.^ ^ The ideal scintillator for gamma ray de- 
tection must have exceptional performance in terms of 
stopping power, luminosity, proportionality, speed, and 
cost. Recently, trivalent lanthanide dopants have re- 
ceived greater attention for fast and bright scintillators. 
In particular, Ce^+ is a favored dopant in many scintil- 
lators due to its allowed optical bd-Af transition which 
is relatively fast (~ 20-40 ns) and it can be doped onto 
La,Y, Gd and Lu sites of many high density host mate- 
rials. Consequently, some of the brightest known scin- 
tillators are Ce-doped such as LaBraiCe,^ LuIaiCe^ and 
Yls'.Ce.^ However, crystal growth an d pr oduction costs 
remain challenging for these materials 

First principles calculations provide a useful insight 
into chemical and electronic properties of materials and 
hence can aid in the search for better materials or guide 
modification of existing materials .1^^ The theoretical 
work presented in this paper is part of a larger project 
^^High-throughput discovery of improved scintillation ma- 
terials^\ which aims to synthesize and characterize new 
materials in microcrystal form and select candidates for 
crystal growth. The main aim of the theoretical studies 
presented here is to develop a fast method to select can- 
didate Ce activated scintillator materials for synthesis as 
well as complement the experimental work through sim- 
ulations of promising synthesized materials. Preliminary 



results from our studies have been presented earlier.1^ 
In this paper we give a detailed account of our first- 
principles calculations and extensive results obtained so 
far using more advanced calculations than presented in 
our previous work. 

The basic mechanism for scintillation in a Ce doped 
material is that an incident gamma ray will produce a 
large number of electron-hole (e-h) pairs in the host ma- 
terial that transfer to the Ce site. The emission of light 
then corresponds to a bd-Af transition on the Ce site 
from the Ce [Xe]4/^5G?^ excited state, usually referred 
to as (Ce^+)*, to the Ce^+ground state [Xe]4f^5d^ (see 
Fig. [T]). Trapping mechanisms on the host, such as self 
trapped excitons, hole traps or electron traps, can quench 
or reduce the transfer of energy to the Ce site (see, for 
example, Ref. [1 for a more detailed discussion of scin- 
tillator mechanisms and quenching processes). 

Recently, there has been a growing interest in ab ini- 
tio calculations of the properties of bd-Af transitions of 
rare-earth ions in solids. Much of this resulted from the 
pioneering work of Dorenbos and collaborators, who in a 
series of papers compiled experimental data of this tran- 
sition and derived semi-empirical models for predicting 
properties of the bd-Af transition and estimated the po- 
sitioning of these states in the host gap.^^ Most of 
the ab initio calculations performed to date for rare- 
earth (RE^+) doping use either cluster models based 
on Hartree-Fock or band structure approaches based 
on Density Functional Theory (DFT). Usually, in the 
embedded cluster calculations, to reduce computational 
costs, the dopant and the first-shell ions around the 
dopant are allowed to relax while the rest of the crystal 
is kept frozen in the crystalline geometry. This can give 
anomalous result&i^'^ and possible deficiencies of this lo- 
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FIG. 1: Schematic diagram of the position of the Ce 
bd and 4/ levels relative to the conduction and valence 
band of the host material for a Ce activated scintillator. 
VBM is the Valence Band Maximum and CBM is the 
Conduction Band Minimum. is the host material 
band gap. 



cal relaxation procedure have been discussed recently by 
Gracia et al.l^for Ce^^ doped YAG (Y3AI5O12). Cluster 
models also cannot give the positions of the conduction 
band (CB) and valence band (VB) of the host relative to 
the dopant states which is closely related to luminescence 
properties. 

In one of the earliest works using a DFT based ap- 
proach, Stephan et al.^^ studied bd-Af transitions for 
a number of trivalent lanthanides using band structure 
calculations. The effect of rare-earth (RE^+) doping in 
semiconducting GaN has also been reported. However, 
within the local density approximation (LDA) or gener- 
alized gradient approximation (GGA) to DFT the self- 
interaction error associated with the localized nature of 
the 4/ shell prohibits the calculation of accurate energy 
differences. There have been attempts to overcome this 
problem using beyond DFT methods, but they have fo- 
cused on studies of bulk Ce compounds. Recently, 
Nishida et al.^^ studied the relationship between the 
local structure around the Ce^+ ion and the emission 
properties of CeFs and Ce203 employing a combination 
of TEM-EELS measurements and first principles band 
structure studies. The 4/-5g? energy gap was shown to be 
in qualitative agreement with known experimental spec- 
tra. 

A band structure approach can also be used to relax 
the doped host matrix to take into account lattice re- 
laxation effects. Andriessen et al.' ^^'^^ ' have performed 
such relaxations for a few known Ce^+ doped scintillat- 
ing compounds. Recently, they have published detailed 
results for the Stokes shift from lattice relaxation studies 
of 4:f-5d excitation of Ce doped lanthanum halide scintil- 
lators using a band structure approach based on density 
functional theory and ionic cluster calculations using the 
Hartree-Fock method.^^ Watanabe et al.l^ studied the 
4/- bd absorption spectra of Ce doped LiYF4 using a 
combination of the pseudopotential plane-wave method 



along with the relativistic molecular orbital approach. 
They found that the Af-bd transitions in the case of Ce^+ 
can be attributed to transitions between molecular or- 
bitals since Ce^+ has a simple [Xe]4/^ electronic struc- 
ture implying that Ce Af-bd transitions can be analyzed 
within the framework of a single-electron approximation. 

Our theoretical calculations for the prediction of can- 
didate scintillator materials are based on studies of the 
Ce 4/ and bd levels relative to the valence band maxi- 
mum (VBM) and conduction band minimum (CBM) of 
the host material, respectively.'!^ A necessary condition 
for scintillation and luminescence is that the Ce 4/ and 
bd levels must be in the gap of the host material. If 
the Ce 4/ level lies in the valence band of the host or 
the bd level is in the conduction band there will be no 
Ce activated scintillation or luminescence. If the 5c? Ce 
state lies below but close to the bottom of the conduction 
band then thermal excitation from the bd state into the 
conduction band can reduce or quench luminescence. It 
should also be noted that under direct optical excitation 
of the 4:f-bd transition some Ce doped systems can show 
strong luminescence but can be weak scintillators due to 
trapping mechanisms on the host that can quench or re- 
duce the transfer of energy from the incident gamma ray 
to the Ce site. 

In the present paper electronic structure calculations 
of Ce- doped compounds are performed with the LDA+U 
(and GGA+U) approach. This method has been shown 
in previous publications to give a better description of the 
localized 4/ states of Ce compared to LDA or GGA.E31I1 
We have tuned the empirical Ueff parameter for the Ce^+ 
impurity atom to match the calculated Ce 4/ to host 
VBM gap with the experimental energy gap for some 
known scintillating and non- scintillating Ce-doped com- 
pounds. Validation and predictions of Ce 4/ -VBM en- 
ergy gaps in the ground state are presented. 

An accurate determination of the Ce bd-CBM en- 
ergy gap for the (Ce^+)* state is difficult using standard 
ground state LDA and GGA approximations to density 
functional theory. A ground state calculation with the 
4/ level filled and the bd level empty yields di bd level 
that will be higher than when the bd level is filled and 
the 4/ level empty. The 4/ level is closer to the nu- 
clei than the bd level so when the 4/ level is emptied 
the screening effect from the positive nuclei will be re- 
duced and the bd level will move lower. The Stokes shift 
can also further lower the bd level but we did not try 
to model that in our simulations. Previous studies have 
found the Stokes shift to be difficult to model accurately 
with DFT based band structure codes. We, therefore, 
performed excited state (constrained LDA) calculations 
and subsequent analysis to allow us to derive a qualita- 
tive measure of the bd-CBM energy gap. Our main aim 
is not an extremely accurate calculation of the bd level 
position, but to determine whether or not it is below the 
CB as this determines if luminescence from the Ce site 
is possible. It should also be noted that the host dopant 
site in our studies is either La, Lu, Gd or Y so the CB 
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has 5c? or 4 d character. Therefore, systematic errors due 
to the LDA type treatment of the Ce 5 c? state wih also be 
present in the determination of the CBM yielding, partic- 
ularly in the case of La,Lu and Gd, a reasonably accurate 
5d-CBM separation due to cancellation of errors. The 
size of the supercells in our calculations typically pro- 
hibited the use of more advanced many- body methods. 
Earlier studies of Ce activated scintillators with the clus- 
ter based Hartree-Fock method found that adding con- 
figuration interaction only had a minor influence on the 
results.'^ Overall we want to develop a high throughput 
method for the screening of large numbers of new ma- 
terials as candidates for bright Ce activated scintillators 
so we restrict our calculations to computationally fast 
first principles methods that can yield good qualitative 
results. 



II. CALCULATION DETAILS 

In order to simulate a dopant in a periodic lattice we 
use the supercell approach with periodic boundary con- 
ditions. We construct a large supercell from periodically 
repeating the unit cell of the host crystal and then re- 
place one of the host trivalent site by a Ce atom. We 
then relax the atomic positions while keeping the cell di- 
mensions fixed. Our basic aim in these studies is to model 
one Ce atom in an infinite host lattice however, the super- 
cell approach introduces spurious dopant-dopant interac- 
tions due to the periodic boundary conditions.^ These 
interactions can cause a broadening of the impurity lev- 
els into bands and also modification of the valence and 
conduction band edges which are the natural reference 
energies for the impurity states. We, therefore, perform 
size scaling studies to be sure the supercells we use are 
large enough to produce converged results for the prop- 
erties of interest. Once we have relaxed the supercell we 
perform a ground state calculation to determine the po- 
sition of the Ce 4/ level relative to the VBM of the host 
material. The filled 4/ level is typically very localized 
and atomic in nature and has almost no bandwidth so 
the 4/-VBM gap is well defined. 

To determine if the (Ce^+)* state lies below the CBM 
we perform a constrained LDA (or GGA) calculation by 
setting the occupancy of the Ce 4/ states to zero and 
filling the first state above the 4/ levels. Previous calcu- 
lations for the (Ce^+)* excited state have been performed 
by removing the Ce 4/ states from the basis functions or 
creating a pseudopotential with Ce 4/ states treated as 
core states. Our method has the advantage that 
we can use the same basis set and pseudopotential for 
both excited state and ground state calculations allow- 
ing direct comparison of energies. We then look at the 
spatial distribution of this excited state to determine if 
it is localized on the Ce or is a delocalized CB charac- 
ter state. The level of localization of an electronic state 
does not have a strict mathematical definition, but for 
the purposes of our studies we will define it as the per- 



centage of the normalized electron density in a Voronoi 
cell centered on the Ce atom. We will also consider rela- 
tive localization: a ratio of localization of a state on the 
cerium site to its next largest localization on a different 
cation (La, Lu, Gd or Y). If the state has no localiza- 
tion (the percentage on the Ce atom is very low and the 
ratio is one or below) then we can consider it is a host 
band structure state and is the bottom of the CB. In 
such a scenario any localized state of Ce 5d character 
lies above the CBM and there is no possibility of scin- 
tillation or luminescence. If the state is localized on the 
Ce and has 5 c? character then we can associate it with 
the so-called (Ce^+)* and a 5c?-4/ transition is possible. 
We found this procedure for determining if there exists 
a (Ce^+)* state below the CBM necessary, as in the sys- 
tems studied there typically seems to be some level of 
hybridization between the host d character CB and the 
Ce 5c? character states. This will, to some degree, delo- 
calize them from an atomic like 5d state centered on the 
Ce. We, therefore, needed a simple way to characterize 
the lowest d type as a CB or Ce state without having to 
resort to very large supercell calculations where the elec- 
tronic states and energies were completely converged. We 
found that for wavefunctions localized over many atomic 
distances the percentage in the Voronoi cell as well as 
the ratio to neighboring cations can be low even though 
to the eye the wavefunction is clearly localized in space. 
Our simple definition of localization does not contain any 
concept of localization distance and poorly characterizes 
states localized over many atomic distances. We will dis- 
cuss this issue more in the results section where we find 
some of the oxide scintillators like YAP (YAlOaiCe) have 
(Ce^+)* states localized over many atomic distances. It 
should also be noted that very localized (Ce^+)* states 
will tend to have a larger binding energy due to Coulomb 
attraction between the 5d electron and the nuclei. The 
removal of the 4/ electron reduces the screening of the 
positive nuclei or can be thought of as leaving a hole state 
(compared to the ground state) which has Coulomb at- 
traction with the filled 5d state. The excited state will 
then have more chance of being lower in energy than na- 
tive exciton states on the host which might otherwise 
reduce or quench the scintillation. The Ce 5d character 
state lying below the CBM is a necessary, but not suffi- 
cient, condition for scintillation. Related to this, one of 
the goals of our work is study how the Ce 5d state prop- 
erties are related to scintillation as well as luminescence 
properties in Ce doped materials. 



A. Atomic relaxation studies 

The initial atomic positions and symmetry information 
of the host crystal were taken from the Inorganic Crystal 
Structure database (ICSD).^^ ^^""The number of atoms in 
the Ce doped supercells was typically 50-150 depending 
on the size of the host unit cell and how many atoms 
were required for reasonable convergence. The Vienna 
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ab initio simulation package (VASPj^^^*^ was used for 
spin-polarized GGA(PBE)^^ and LDA calculations. The 
projector-augmented wavefunction (PAW) approach, de- 
veloped by Blochl^^ and adapted and implemented in 
was used for the description of the electronic 
wavefunctions. Plane waves have been included up to an 
energy cut-off of 500 eV. Integration within the Brillouin 
zone was performed with a F point centered grid of k- 
points. The number of irreducible A;-points was typically 
chosen to be 4 or 8 depending on the size and geometry 
of the supercell. The energy convergence criterion was 
set to 10~^ eV and the maximum component of force 
acting on any atom for relaxation of the atomic positions 
after doping with cerium was checked to be less than 
0.01 eV/A in every direction. Cerium pseudopotential 
was chosen to include (55,5p,65,4/,5(i) as valence elec- 
trons. We have used the rotationally invariant method of 
Dudare\^^ as implemented in VASP^^ for an on-site +U 
correction to treat the cerium 4/ electrons with a single 
parameter Ueff = U-J, where the Hubbard U parame- 
ter is the spherically averaged screened Coulomb repul- 
sion energy required for adding an extra electron to the 
Ce 4/-states and the parameter J adjusts the strength 
of the exchange interaction. We determined Ueff em- 
pirically by adjusting it to correspond to experimental 
results (see Section IV. A). An artifact of DFT-PBE (or 
LDA) calculations is that unoccupied La 4/ states are po- 
sitioned at the bottom of conduction band.l^ However, 
La 4/ states lie higher in energ}'^^*^ so in our calcula- 
tions we push the La 4/ states higher in the energy plot 
using the the LDA-hU approach with the Ueff parameter 
taken from Ref. [46]. Calculations for Gd systems were 
performed with the Gd_3 (4/ states in the core) pseu- 
dopotential. We checked two test calculations with the 
regular Gd pseudopotential (4/ electrons as valence) and 
found the results to be very similar. We do not employ 
spin-orbit coupling in our calculations as in the case of 
the La halides this was found to only move the Ce 5d 
states by a maximum of about 0.2 eW^ which would not 
change our qualitative conclusions and would increase the 
computational cost. 

For the purposes of comparison and checking the ac- 
curacy of the pseudopotentials, ground state density of 
states (DOS) calculations were done for a few Ce-doped 
systems using the full potential linear augmented plane 
wave (FP-LAPW) code WIEN2K.'^The relaxed atomic 
positions from the VASP code were used as input to 
the WIEN2K code. The same GGA(PBE) functional 
was used in the two codes. We kept the A:-point grid 
and energy convergence criteria similar to VASP calcu- 
lations. The number of plane waves was restricted to 
Rmt X kmax = 7. The fully localized limit (FLL) form of 
GGA+U implementation was used within the WIEN2K 
code to treat the Ce 4/ orbital. The value of the Ueff pa- 
rameter was kept the same in the two calculations. The 
results were found to be very similar to PAW calcula- 
tions with VASP with the positions of the various bands 
varying by only a few percent between the two codes. 



B. Excited state calculations 

In order to study the (Ce^+)* state constrained LDA 
(and GGA) calculations were done at the F point using 
the VASP code. The occupation numbers were manually 
set to empty the Ce 4/ states and fill the next highest 
state. The band decomposed charge density was subse- 
quently analyzed to derive the localization parameters. 

Excited state calculations were also done within 
the PAW framework as implemented in the ABINIT 
code.'^^H^ The electronic wave functions were expanded 
in plane waves up to a kinetic energy cutoff of 60 Hartree. 
Self-consistency was achieved using a A:-point grid cen- 
tered at the F point in reciprocal space. The energy tol- 
erance for the charge self-consistency convergence was set 
to 1 X 10~^ Hartree. Band-decomposed charge density at 
the lowest energy A;-point was subsequently analyzed to 
derive the localization parameter. ABINIT calculations 
were, however, limited to compounds with elements hav- 
ing reliable PAW data sets. 

III. THEORETICAL CRITERIA FOR 
SCINTILLATION AND LUMINESCENCE 

Based on present understanding of scintillation physics 
and our previous first-principles studies of known Ce^+ 
scintillators (e.g., YAlOsiCe (YAP), LusSiOsiCe (LSO), 
LaCl3:Ce, LaBraiCe, LalgiCe, LusSisOriCe (LPS) etc.) 
and non Ce- activated scintillators (e.g., Y203:Ce, 
La203:Ce, LaAlOaiCe (LAP)) we have developed three 
criteria based on the following theoretically calculable pa- 
rameters to predict candidate materials for bright Ce^+ 
activated scintillation.'^ 

1. The size of the host material bandgap. 

2. The energy difference between the VBM of the host 
and the Ce 4/ level. 

3. The level of localization of the lowest d character 
excited state needed to determine if it is a host CB 
state or a Ce 5(i character state. 

Criterion 1 is related to the fact that the number of 
electron- hole pairs produced by an incident gamma ray 
is inversely proportional to the bandgap energy although 
the constant of proportionality varies from material to 
material.!^ Therefore, the band gap should be as small 
as possible but must be large enough to accommodate 
the Ce 4/ and 5d states. LDA and GGA are known to 
underestimate the bandgap, but for the purposes of our 
calculations it does provide trends in families as well as 
comparative results for similar materials. More accurate 
bandgaps can be calculated theoretically by using more 
advanced methods that go beyond LDA but these meth- 
ods are typically more computationally costly. Therefore 
for the purposes of a qualitative prediction of candidate 
scintillator materials we use LDA and GGA calculations 
of bandgaps. 
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Criterion 2 is related to the cases where the energy 
transfer to the Ce site occurs by sequential hole trapping 
and electron trapping on the Ce site. For these cases if 
the 4/-VBM gap is large there will be a low probability 
of the hole transferring from the host to the Ce site via 
thermal excitation which will reduce scintillation bright- 
ness. 

Criterion 3 as discussed in the previous section is how 
we determine if the lowest d character excited state can 
be associated with a Ce 5g? character state or a con- 
duction band state of the host material depending on 
whether the state is localized or not. 

A further expected result of our calculation may be 
estimation of the Ce bd-GBM. gap. This is in fact ex- 
tremely difficult to calculate accurately as in our large 
supercells there are many d character bands associated 
with the host as well as the those associated with Ce. 
The bands associated with Ce localized bd states, even 
with large supercells, typically still have some curvature 
due to finite size effects making the 5(i-CBM gap not well 
defined. It is also a difficult process to scan up through 
the lowest d character bands to determine which are as- 
sociated with Ce bd states and which are CB character 
states as most bands show some level of hybridization 
particularly at higher energies and presumably close to 
to the CBM. In many cases the excited state calculation 
is also slow or problematic to converge due to the close 
proximity of the many d character bands to the filled 
d band. This is also one of the reasons we have used 
different codes such as ABINIT and VASP for these cal- 
culations as we have found that for different systems one 
code may have better convergence properties than the 
other due to the different minimization methods used. 
For some known scintillators and non-scintillators we 
have performed more detailed studies of the character 
of the different d bands and we will present this in future 
work. In particular for some non-scintillators we do find 
bd character Ce states within the conduction band al- 
though they typically have some hybridization with the 
host d states. The problems in calculating a 5(i-CBM 
gap are not shared in determining the 4/-VBM gap as 
the 4/ state is typically extremely localized on the Ce 
atom giving a flat band even with modest sized super- 
cells. The higher energy empty 4/ states are also well 
separated from the filled 4/ states giving fast convergence 
for the ground state calculation. 

From the point of view of predicting new bright scin- 
tillators it would be useful to develop more theoretical 
criteria related to trapping processes on the host that 
can limit energy transfer to the Ce site. Unfortunately, 
it is difficult to use first principles calculations to develop 
criteria related to these host processes as the exact na- 
ture of the trapping sites on the host are often poorly 
understood from experiment as well as the details of the 
energy transfer mechanisms to the Ce site. For example, 
accurate calculations of deep self-trapped host excitons 
such as those found in LaFa often require advanced many 
body theories and involve significant lattice relaxation. 



The dynamical nature of the transfer processes of host 
excitons and hole or electron traps to the Ce site is also 
difficult to model from first principles although there has 
been work done in developing empirical models of these 
processes. Overall though, from an energetics point of 
view, we would expect the transfer of energy to the Ce 
site to be most favored, the deeper the (Ce^+)* state is 
within the bandgap of the host material. For example, if 
the (Ce^+)* state is lower in energy than any host STEs, 
it will preferentially form provided there are no large en- 
ergy barriers to transfer processes from the host to the Ce 
site. As with semiconductor dopant states the depth of 
the (Ce^+)* state in the gap of the host material will be 
related to its level of localization.!^ This is particularly 
true where the character of the CB and dopant state are 
the same which is the case for the systems studied here. 
Hence we expect criterion 3 to also be related to the 
brightness of a Ce activated scintillator. 



IV. RESULTS AND DISCUSSIONS 

In this section we present results of our theoretical 
studies on Ce doped compounds. The discussion is di- 
vided into three subsections. The first subsection con- 
cerns the ground state density of states calculations for 
Ce^+ doped compounds specifically, the determination 
of the Ueff parameter from experimentally measured Ce 
4/-VBM energy gaps. Cell size scaling studies were also 
performed to check the dependence of the Ueff parameter 
on the simulation supercell size. The second subsection 
presents results of the excited state calculations. Simu- 
lation cell size scaling studies are also presented in this 
section. In the last subsection we perform calculations for 
some new materials doped with Ce and apply our theo- 
retical criteria for the prediction of new bright candidate 
Ce activated scintillators. One of the new scintillators 
predicted by our calculations was synthesized in micro- 
crystal form and confirmed to be relatively bright. We 
have also generated a database of Ce 4/-VBM energies 
predicted from first-principles calculations of more than 
100 compounds. 



A. Ground state calculations : Determination of 
Ueff parameter 

Ueff can be determined in a self-consistent way as 
demonstrated by Cococcioni et al.^^ for Ce02. However, 
more frequently Ueff is chosen in such a way as to repro- 
duce with reasonable accuracy an experimentall y pea- 
sured quantity like cell volume, bulk modulus, etc.'^^'^^IIl 
We are not aware of any prior publication related to 
LDA+U type calculations for Ce-doped insulators of the 
type used for scintillator detectors. Fortunately, experi- 
mental measurements of Ce 4/-VBM gap are known for a 
few scintillators so we chose the Ueff parameter to closely 
match these known gaps. 
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Figure |2] shows the total density of states plot for 
LaBrsiCe for different values of the Ueff parameter from 
a GGA(PBE)+U calculation. The filled Ce 4/ state is 
at the Fermi level which is set to zero. The experimen- 
tally measured Ce 4/-VBM gap for LaBraiCe is 0.9 eV 
(±0.4 eV).^^ We observe from the figure that the calcu- 
lated 4/-VBM energy gap using Ueff = 2.5 eV matches 
the experimental data. It is important to note that the 
value of Ueff used in the literature for bulk Ce(III) com- 
pounds (Ueff = 4.5 eV for PBE functional^"^ is differ- 
ent from our results. This is mainly due to the itinerant 
nature of 4/ electrons in Ce bulk compounds which par- 
ticipate in bonding compared to our doped ionic systems 
where the single Ce 4/ electron is atomic-like in nature. 
Hence it is important to tune the empirical parameter 
Ueff to get a close match with experimental data for Ce- 
doped scintillator materials. It should also be noted in 
these DOS plots that the filled Ce 4/ state is close to 
being a delta function corresponding to a fiat band. The 
4/-VBM gap is therefore well defined from band struc- 
ture calculations for these types of system. 
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FIG. 2: Ground-state DOS plot for LaBraiCe from 
PBE+U spin polarized calculations for different values 
of Ueff. Fermi level is set at 0. U-J =0 eV corresponds 
to DFT-PBE result. Experimentally estimated 
4/-VBM gap is 0.9 ± 0.4 eVW^ 

Figure |3] shows total DOS plots of LusSisOriCe (LPS) 
for different simulation cell sizes for a fixed Ueff param- 
eter. Even for these small cell sizes there is negligible 
variation in the Ce 4/-VBM energy gap (~2%) with cell 
size. All the data presented in Table |T] is for similar or 
larger cell sizes so we are confident that any finite size 
effects on the Ce 4/-VBM energy gap are below a few 
percent. However, as we show in the next subsection, 
simulation cell size has a greater infiuence on the local- 
ization of the excited state. 

We repeated the calculations for a few systems such 
as YlsiCe and LaBrsiCe using the LDA+U functional 




FIG. 3: Density of states (DOS) plot for Lu2Si207:Ce 
from PBE+U (U-J = 2.5 eV) spin polarized 
calculations for two different cell sizes. 



and found negligible change in the results compared to 
GGA(PBE)+U calculations. The choice of the approx- 
imation to the exchange-correlation functional (PBE or 
LDA) does not affect the position of the impurity (Ce) 4/ 
levels. This is unlike calculations for bulk Ce compounds 
where different values of Ueff have been used for differ- 
ent functionals.'^ This is because the Ce 4/ atomic-like 
character changes little for the different functionals. 

Table [J summarizes the results of our studies to 
tune the Ce 4/-VBM to known experimental measure- 
ments for Ce doped materials. We see from Table [T| 
that GGA(PBE)+U calculations with a Ueff = 2.5eV 
give good agreement with experimentally measured Ce 
4/-VBM gaps for most materials with the exception of 
LalaiCe where Ueff = 2.2 eV gave the best agreement 
with experiment. LalaiCe is one of the smallest bandgap 
scintillator materials so the bonding is more covalent in 
nature than in other scintillator materials. This may ac- 
count for the slightly different character of the Ce 4/ state 
requiring a 0.3 eV lower value of Ueff than in the other 
systems. For many of the experimental results reported 
in the table, error bars are not quoted in the publications. 
We have found that in scintillator materials the charac- 
ter of the Ce 4/ is extremely atomic and very similar for 
different hosts which explains the universality of the Ueff 
value in this class of materials. In the case of the heavier 
host La halides there may be some weak dependence of 
the character of the Ce 4/ on the local environment sur- 
rounding it. From Table [l| we can see that the 4/-VBM 
gap is larger for fiuorides and progressively decreases as 
we go down the periodic table from F to I. Oxides typi- 
cally have a larger 4/-VBM gap than sulfides. The Ueff 
parameter we found that gives the best fit to experiment 
has no variation from oxides to halides and thus provides 
a potentially simple method to find the Ce 4/-VBM gaps 
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TABLE I: Experimentally measured and calculated 
(PBE+U) 4/-VBM gaps for known Ce-activated 
scintillators and non-scintillators. 



Compound 




± l3 12j-\- \J 




gap (eV) 


result (eV j 


iJcxSDL 3 . v_yc 


n Q _|_ n 4 

U.c/ ZL U.^ 


n Q 

u.c/ 


(scintillator) 


(Dorenbos et al.^ 


(Ueff = 2.5 eV) 


T ill O S1 1 O O '7 ■ (jP 


2.9 


2.9 


(scintillator) 


(Pidol et alP) 


(Ueff = 2.5 eV) 


T.iioSiOt^-Cp 


3.1 


2.7, 2.9 


(scintillator) 


(Joubert et al.^ 


(Ueff = 2.5 eV) 






two substitu- 






tion sites 


YAlOsiCe 


- 3.3 


3.0 


(scintillator) 


/at:i 1 -^1 160|\ 

[NiKL et al.' — 'j 


(Ueff = 2.5 eV) 


LalsiCe 


0.2-0.3 


0.25 


(weak scmtiiiatorj 


/T> • ^1 I61K 

(Bessiere et al.' — ■) 


(Ueff = 2.2 eV) 


YP04:Ce 


- 4.0 


3.65 


(weak scintillator) 


/TA „ 1 J62|\ 

(JJorenbos — 'j 


(Ueff — ^.0 eV ) 


LaAlOsiCe 


- 2.0 


2.1 


(non-scintillator) 


(van der Kolk et al.^^ 


(Ueff = 2.5 eV) 


La203:Ce 


- 2.8 


2.9 


(non-scintillator) 


(Yen et al.^^ 


(Ueff = 2.5 eV) 


YsOsiCe 


- 3.4 


3.4 


(non-scintillator) 


(Pedrini et al.'^ 


(Ueff = 2.5 eV) 



for different types of compounds as compared to precise 
measurements.^^ In our calculations reported in subse- 
quent sections for new materials we used Ueff = 2.5 eV 
to correct the 4/ position except for iodides and sulphides 
where we used 2.2 eV. 

In all the systems studied we found the 4/ level to 
be above the VBM so unlike the 5d level relative to the 
CBM, the 4/ level position relative to the VBM is not a 
factor in preventing 5d-4:f emission in Ce doped systems. 

Our studies also revealed that ionic relaxation was pre- 
dominantly influenced by the difference between Ce^+ 
ionic radii and the trivalent host cation dopant site with 
the choice of Ueff parameter having negligible effect. Ce- 
doped Lu^+ and Y^+ compounds showed significant re- 
laxation as compared to La^+ compounds primarily be- 
cause of the almost 10% size mismatch between Ce^+ and 
Lu^+ , Y^+ and less than 1% mismatch between the ionic 
radii of Ce^+ and La^+ . 



B. Excited state calculations 

As described in Section |IIB[ we performed excited 
state calculations by manually setting the occupation of 
all the Ce 4/ states to zero and filling the next highest 
state. Figure |4ja-b) shows the atom projected partial 
density of states for Lu2Si207:Ce (LPS) in the ground 
state and excited state. There is no atomic relaxation in 
the excited state so there is no Stokes shift in our cal- 
culations. For this system the valence band of the host 
material consists of O p states hybridized with Lu 4/ 
states and the conduction band consists of Lu 5d char- 
acter states. We can see from the ground state plot that 
there are Ce 5d states below the Lnbd states even in the 
ground states DOS and these move about 0.5 eV lower 
relative to the CBM in the excited state plot. In the 
excited state plot the excited-fermi level lies above the 
lowest Ce 5d level showing the filling of the lowest Ce bd 
level. In this system the lowest Ce 5d levels are clearly 
below the host CB as is necessary for the Ce bd-Af tran- 
sition to occur. As can be seen from this plot the Ce 
5d levels have some bandwidth resulting in a continuous 
DOS function for the Ce 5d states rather than the delta 
type function we find for the very localized Ce 4/ states. 
The higher energy Ce 5 c? character states are hybridized 
with the Lu 5c? states. 

Figure [5] shows the atom projected partial density of 
states for Ce doped LaBra in the ground state which, un- 
like LPS, is more typical of the type of result we obtained 
for different scintillators. The 5d states on the La and Ce 
are hybridized and occur at the same energy so there are 
no well defined Ce character 5d states below the CB. For 
these types of systems we find the characterization of the 
lowest filled excited d state in terms of its localization on 
Ce to be the best method to determine if it has Ce 5 c? 
character or is a host CB character state. 

Figure [6] shows charge density isosurface plots of the 
first d character excited state at the gamma point 
for some known scintillators and non-scintillators. For 
the known non-scintillators La203:Ce, Y203:Ce and 
LaAlOaiCe there is no localization on the Ce and the 
excited state has a band structure character distributed 
throughout the supercell. On the other hand, a local- 
ized excited state with d character forms on the Ce 
site for the known scintillators Lu2Si207:Ce, LaBraiCe 
and YAlOsiCe. As can be seen from the plots there is 
a large range of localization of the excited state with 
Lu2Si207:Ce being much more localized than the other 
systems. 

Table [ll] presents a list of our theoretically calculated 
parameters of bandgap, 4/-VBM gap, % localization and 
localization ratio for a list of known scintillators and non- 
scintillators compared to experimental data for bandgaps 
and scintillation luminosity. As expected LDA consis- 
tently underestimates the bandgap but does correctly 
predict the ordering of bandgaps for similar materials 
and families of materials. In all these materials the Ce 4/ 
level is above the VBM so the occurrence of the 4/ level 
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FIG. 4: Atom projected partial density of states (DOS) 
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in the ground state (a) and excited state |(b)] Fermi 



level is set to 0. / character states are shown in red, d 
states in blue, p characters states in green and s 
character states are shown in black. 



within the VB never seems to be a factor in quenching lu- 
minescence and scintillation in Ce doped materials. Also, 
to the best of our knowledge, there is no experimental ev- 
idence of Ce 4/ states inside the host VB. The main result 
from this table is that we have essentially no localization 
of the lowest excited d state for all the non- scintillators. 
The brightest scintillators typically have low bandgaps 
and small 4/-VBM gaps although it should be noted that 
for scintillation the bandgap has to be large enough to 
accommodate the Ce 4/ and 5d states. Overall, there 
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FIG. 5: Atom projected partial density of states plots 
for LaBraiCe in the ground state. Fermi level is set at 
eV. Calculation used the (GGA)PBE functional and 
Abinit code 



is good qualitative agreement between our three criteria 
and bright scintillators. 

The La halides represent a family of materials that 
have been very heavily studied experimentally for Ce ac- 
tivation as they are all scintillators and have a large range 
of bandgaps. LalaiCe has a very low band gap of 3.3 eV 
and is thermally quenched at room temperature due to 
the proximity of the excited Ce state to the CBM but 
has reasonable luminosity at lOOK.^^ Excited state cal- 
culations for this system are particularly difficult to con- 
verge since Ce bd states hybridize and are very close to 
the host CB. This is consistent with the experimentally 
estimated bd-CBM gap of - 0.2 eV.I^ This also leads 
to a relatively low values for the % localization and ra- 
tio. The (Ce^+)* excited state is favorably localized for 
LaClsiCe and LaBraiCe. This agrees with the fact that 
these materials are well known bright scintillators used 
in several gamma ray detection applications.!^ LaBr 3: Ce 
in particular has lower bandgap, favorable 4/-VBM gap 
and reasonable localization on the Ce site. It should be 
noted that the role of host STEs is known to be impor- 
tant in the transfer of energy to the Ce site for LaClaiCe 
and LaBrsiCe where the transfer mechanism is efficient 
leading to bright scintillators. In these cases the size of 
the 4/-VBM gap will play less of a role in determining 
the brightness .I^LaFaiCe is an example of a system that 
is known to have a very deep STE of a lower energy than 
the (Ce^+)* excited state.^^ Even though the lowest d 
character excited state is of Ce 5d character and well lo- 
calized this limits the transfer of energy to the Ce site 
and results in very low luminosity for this particular ma- 
terial. Moreover, the bandgap and 4/-VBM energy gap 
for this system are quite large which leads to compara- 
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FIG. 6: Lowest character excited state plots for Ce scintillators and non-scintillators. Plots show charge density 
isosurfaces of the excited states . Ce atom is shown in blu e, ra re-eart h io n (=La, Lu, Y) is in yellow and the anions 
are shown in red. [(a)] LaBrs ; [(b)] Lu2 Si2 O7 ; [(c)] YAIO3 ; p)] Y2 O3 ; [(e)] La2 O3 ; [(f)] LaAlOg . The excited state is 



delocalized (very little or no concentration around Ce site) for non-scintillating compounds La203:Ce, Y2 03:Ce and 
LaA103:Ce. However, Ce^+ scintillators have good localization of the excited state on the Ce site. 



tively lower e-h pair production and a low probability of 
sequential hole and electron capture by Ce. Thus, even if 
there were no low energy host STEs we would not expect 
this system to be a bright scintillator. 

Oxide scintillators in general have wider bandgaps so 
Ce 4/ and 5d states are mostly better separated from the 
band edges. As we can see from Table [n] the {Ce^^y 



state is favorably localized in most of these systems. 
YA103:Ce is an example of a system that has a rather 
low % localization and ratio even though Figure [g] clearly 
shows a localized state. The main reason for this is that 
the state is localized over a few interatomic spacings so 
since our simple measures of localization have no mea- 
sure of localization with distance from the Ce they tend 
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to give low values for these types of localized states. 

We also studied the dependence of localization of the 
excited state with cell size. As an example, Figure [7| 
shows the localization of a bright scintillator YlsiCe with 
increasing cell size. YI3 has a trigonal crystal structure 
with a= 6=7.4864 A and c=20.88 A. In the excited state 
plot for the 24 atom conventional unit cell, even though 
the dopant Ce^+ ion has a high percentage localization, 
the ratio to the next highest indicates that Ce sites 
in the periodically repeated cells are interacting with each 
other in the direction of the shortest cell dimension (hori- 
zontal plane containing the Ce atom and the neighboring 
Y). Now when we scale in the horizontal dimensions for 
the 96 atom simulation cell we find that the Y atoms 
in the same plane as Ce have some fraction of the ex- 
cited state, but Ce has the highest percentage of the 
localization of the excited state. This, still, does not 
clearly show a predominating Ce localization expected of 
a bright scintillator like YlaiCe because the excited state 
wavefunction is not well localized within the cell volume 
and consequently, there is interaction with the Ce sites 
in the periodically repeated cells in the plane contain- 
ing lattice vectors a and b. Upon scaling the simulation 
cell size to 384 atoms the excited state becomes predom- 
inately concentrated on the Ce site. We found that the 
convergence with cell size varied significantly for differ- 
ent host materials. The cell sizes quoted in Table [ll| were 
chosen to give well converged results for the materials 
studied and are typically smaller than for YlaiCe. 



C. Prediction of new candidate Ce scintillators 

The next phase of our theoretical studies was to ap- 
ply the criteria we have developed from studying known 
scintillators and non-scintillators to the prediction of new 
candidate scintillators. In this study we performed the 
same calculations for about a hundred new materials as 
we did for the known materials. We have chosen new 
host compounds based on their stopping power and their 
ease of doping with Ce (i.e., the availability of a trivalent 
sites such as Y, La, Lu or Gd for substitution by a Ce 
atom). In particular, we studied the BaY^X^iCe (X= 
F, CI, Br, I) family of materials where BaY2F8:Ce is known 
to be a weak scintillator but the performance of the other 
materials doped with Ce was unknown. The new mate- 
rials we have studied having the best characteristics for 



bright Ce activation are listed Table III 



We found very strong localization of the (Ce^+)* state 
for Ba2YCl7:Ce (see Fig.[8|. The bandgap and 4/-VBM 
separation for this material have values that are close 
to those of some of the well known scintillators such as 
LaClsiCe. On the basis of our theoretical criteria out- 
lined in Section |III[ Ba2YCl7:Ce was expected to be a 
good candidate for a bright new scintillator. It was subse- 
quently synthesized and found to be bright in microcrys- 
tal form.^ In terms of predictions of non-scintillators we 
have studied many other families of materials and have 



found that for all Y and La host materials containing Ti, 
Zr and Hf there is no localized excited Ce state below the 
conduction band. Ce^+ doping in Bi^+ host compounds 
also leads to no localized Ce 5d state below the conduc- 
tion band. Some of these studies will be the subject of 
future publications. We have also previously published 
theoretical work on Ce doped Y and La oxyhalides^^ as 
well as Y halides^^ which included known as well as new 
scintillators. 



V. CONCLUSIONS 

In this paper we have presented DFT based first prin- 
ciples studies for Ce activated scintillator detectors. The 
main aim of this work was to determine what theoreti- 
cally calculable parameters are easily related to lumines- 
cence and scintillation. To more accurately calculate the 
4/-VBM position we used the LDA+U approach where 
we determined Ueff by comparison with experimental re- 
sults. We found that a value of Ugff of 2.5 eV gave good 
agreement with experiment for a wide range of scintil- 
lator materials. Based on this we have calculated the 
4/-VBM gap for many known and new materials, some 



of which are presented in Tables [TT| and III We have also 
generated a database of Ce 4f -VBM energies predicted 
from first-principles calculations for more than a hundred 
new compounds. We also performed excited state calcu- 
lations using a constrained LDA approach to determine if 
the first excited d character state was localized on the Ce 
or was of conduction band character. From these stud- 
ies we developed a set of theoretically calculable criteria 
that characterize bright Ce scintillation. We then vali- 
dated these criteria by studying known scintillators and 
non-scintillators. These criteria were then calculated for 
about a hundred new materials to determine if they were 
candidates for bright Ce activation. The best candidates 



are listed in Table III This approach, involving first- 
principles calculations of modest computing requirements 
was designed as a systematic, high- throughput method to 
aid the discovery of new bright Ce activated scintillator 
materials. This approach has also been extended to Eu 
and Pr doped systems which will be reported in future 
publications. 
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TABLE II: Calculated DFT-PBE bandgaps and energy differences for known Ce activated scintillators and 
non-scintillators. Experimental luminosity data in photons/MeV is taken from Ref. [67 and the references therein. 

** corresponds to no observed Ce emission. 



Compound 

(atoms in supercell) 


LDA Bandgap^ 

(eV) 


Ce 4/-VBM gap 

(eV) 


(Ce3+ 

fry 


)* localization 
ratio 


Luminosity 

(photons/Me V) 


LaFs (48) 


7.8 (9.7^^) 


3.5 


46 


9.14 


2200 


LaCls (128) 


4.6 (7*^ 


1.4 


40 


6.08 


48000 


LaBrs (128) 


3.6 (5.9^^) 


0.9 


21 


5.70 


74000 


Lais (64) 


1.6 (3.3^^) 


0.25 


18 


2.52 


200-300^ 


LaMgBsOio (68) 


5.7 (8.8) 


2.6 


18 


2.48 


1300 


YI3 (384) 


2.8 (- 4.1^ 


0.6 


31 


3.48 


98600 


YAIO3 (160) 


5.4 (8.5-8.4^ 


3.0 


21 


3.17 


21600 


LiGdCU (96) 


4.6 


1.4 


74 


27.6 


64,600 


LusSisOr (88) 


5.5 (7.8^^) 


2.9 


55 


6.8 


26000 


LusSiOs (64) 


4.8 (6.6^^) 


2.9 


33 


7.3 


33000 


CssLiYCle (40) 


5.0 (> 5.d^ 


1.8 


50 


5.8 


21600 


/3-KYP2O7 (88) 


5.9 (-7.71^ 


2.7 


35 


6.4 


10000 


LaAlOs (120) 


4.0 (5.^ 


2.1 


4 


1.6 


** 


Y2O3 (80) 


4.6 (5.#^ 


3.4 


2 


0.63 


** 


LasOs (40) 


4.0 (5.3-5.#^ 


2.9 


1 


0.15 


** 


LU2O3 (80) 


4.7 (5.8) 


2.9 


2 


1.1 


** 


GdsOs (80) 


4.4 (5.4) 


2.8 


4 


0.9 





^The value in parentheses refers to known experimental bandgaps. 
^Luminosity 16000 ph/MeV at 100 K. 
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FIG. 7: Excited state charge density isosurface plots of Yl3:Ce showing the effect of scahng the simulation cell size. 
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TABLE III: Calculated bandgaps, 4/-VBM separation and localization for new Ce doped compounds. 
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